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ABSTRACT 

To evaluate the effect of turbulent heating in the thermal balance of interstellar clouds, we develop 
an extension of the log-Poisson intermittency model to supersonic turbulence. The model depends on 
a parameter, d, interpreted as the dimension of the most dissipative structures. By comparing the 
model with the probability distribution of the turbulent dissipation rate in a simulation of supersonic 
and super-Alfvenic turbulence, we find a best-fit value of d = 1.64. We apply this intermittency model 
to the computation of the mass- weighted probability distribution of the gas temperature of molecular 
clouds, high-mass star-forming cores, and cold diffuse HI clouds. Our main results are: i) The mean 
gas temperature in molecular clouds can be explained as the effect of turbulent heating alone, while 
cosmic ray heating may dominate only in regions where the turbulent heating is low; ii) The mean gas 
temperature in high-mass star-forming cores with typical FWHM of ~ 6 km s^ 1 (corresponding to a 
ID rms velocity of 2.5 km s _1 ) may be completely controlled by turbulent heating, which predicts a 
mean value of approximately 36 K, two to three times larger than the mean gas temperature in the 
absence of turbulent heating; iii) The intermittency of the turbulent heating can generate enough hot 
regions in cold diffuse HI clouds to explain the observed CH + abundance, if the rms velocity on a scale 
of 1 pc is at least 3 km s _1 , in agreement with previous results based on incompressible turbulence. 
Because of its importance in the thermal balance of molecular clouds and high-mass star-forming 
cores, the process of turbulent heating may be central in setting the characteristic stellar mass and in 
regulating molecular chemical reactions. 

Subject headings: ISM: kinematic and dynamics — turbulence 
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1. INTRODUCTION 

The dissipation of turbulent kinetic energy provides 
a potentially important heating source in a variety of 
Galactic astrophysical environments, such as the solar 
wind (e.g., Matthaeus et al. 1999), interstellar clouds 
(e.g., Falgarone and Puget 1995), and the warm ionized 
medium (e.g., Minter and Balser 1997). The effect of 
turbulent heating has also been studied in extragalactic 
environments, such as in the context of the broad-line 
regions of quasars (Bottorff & Ferland 2002), and in in- 
tracluster cooling flows (Dennis & Chandran 2005). 

Astrophysical turbulence is often highly supersonic and 
magnetized. Energy decay in supersonic turbulence was 
thought to be very rapid due to shocks (e.g., Goldre- 
ich and Kwan 1974), while the presence of strong mag- 
netic fields was believed to suppress the compressible 
modes and thus reduce the dissipation rate (Arons & 
Max 1975). Recent numerical simulations offered a bet- 
ter understanding of energy dissipation in supersonic 
MHD turbulence (Stone et al. 1998, Mac Low et al. 
1998, Padoan and Norlund 1999). As in incompressible 
hydrodynamic turbulence, where the kinetic energy de- 
cays in a turnover time of the largest eddies, the dissipa- 
tion timescale in supersonic MHD turbulence is of order 
the flow crossing timescale (Stone et al. 1998) or, equiva- 
lently, the dynamical timescale at the driving scale (Mac 
Low 1999). 

In typical molecular clouds the dynamical timescale is 
rather short, of order a million year, suggesting the need 
of continuous energy injection to support the observed 
turbulence. This result also implies a large turbulent 
heating rate. An estimate of the average dissipation rate 



from this timescale shows that it can be several times 
larger than the cosmic-ray heating rate, and thus may be 
the primary heating source in molecular clouds. Stone et 
al. (1998) also argued that the average turbulent heating 
rate can be comparable to the photoelectric heating in HI 
clouds with large velocity dispersions (however, accord- 
ing to Wolfire et al. 2003, the overall turbulent heating 
in the neutral medium may not be sufficient to produce 
the observed CII luminosity). In this paper, we provide 
a general theoretical formulation of the problem of tur- 
bulent heating, and investigate its effect on various types 
of interstellar clouds. 

Turbulent dissipation is characterized by its strong 
spatial roughness (see Figures [1] and [5]) . Extreme dis- 
sipation events appear in the smallest structures occu- 
pying a tiny volume or mass fraction, while a significant 
fraction of the flow experiences essentially no dissipation. 
This implies a broad probability distribution of the dis- 
sipation rate, which must be taken into account for a 
consistent investigation of turbulent heating. The ex- 
tended tail of this distribution at large dissipation rates, 
corresponding to highly dissipative structures, is respon- 
sible for the anomalous scaling of the high-order velocity 
structure functions, referred to as intermittency in tur- 
bulence theory (Frisch 1995). 

Intermittency has been extensively studied in incom- 
pressible turbulence. The intermittent model by She and 
Leveque (1994), which considers a hierarchy of dissipa- 
tion rates of different levels and relates them to the frac- 
tal dimension of the most intermittent dissipative struc- 
tures, has been very successful in reproducing the mea- 
sured scaling exponents of structure functions in incom- 
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pressible turbulence. It has been shown that the model 
is equivalent to a log-Poisson distribution of the dissipa- 
tion rate (Dubrulle 1994, She and Waymire 1995). The 
probability distribution of the dissipation rate in super- 
sonic turbulence has not been studied yet. Although the 
She and Leveque model with a fractal dimension of 2, 
corresponding to shocks, agrees well with the structure 
functions computed from numerical simulations of highly 
compressible turbulence (Boldyrev et al. 2002, Padoan 
et al. 2004), it remains to be confirmed whether the dis- 
tribution of the dissipation rate in supersonic turbulence 
is consistent with a log-Poisson process. This theoreti- 
cal concern and the wide application of an intermittency 
model of turbulent heating to various astrophysical en- 
vironments are the primary motivations of the present 
work. 

An important effect of the intermittency of the turbu- 
lent dissipation is the generation of small regions with 
very large heating rate, and thus very high temperature. 
Falgarone and Puget (1995) were probably the first to 
recognize the importance of this effect in cold HI clouds. 
Observed molecules, such as CH + , suggest the existence 
of hot regions in cold HI clouds, because their produc- 
tion requires temperatures much higher than the aver- 
age. Adopting experimental results from incompressible 
turbulent flows, they found that strong local turbulent 
heating in HI clouds could produce a sufficient fraction 
of hot regions to explain the observed abundance of CH + 
molecules. In this work we address the validity of their 
result in the case of an intermittency model for super- 
sonic turbulence, more appropriate for cold HI clouds. 

In §2, we study the intermittent energy dissipation in 
supersonic turbulence. We show that the log-Poisson 
intermittency model gives probability distributions for 
the dissipation rate in excellent agreement with those 
at the resolved scales of numerical simulations of super- 
sonic turbulence. In §3, we discuss heating and cooling 
processes in the interstellar medium and give the energy 
balance equations. We apply the log-Poisson intermit- 
tency model to investigate turbulent heating in molec- 
ular clouds, high-mass star-forming cores, and cold HI 
clouds in §4. Discussions and conclusions are given in 
§5- 

2. INTERMITTENCY OF TURBULENT HEATING 

The energy dissipation rate in turbulent flows is known 
to exhibit fluctuations in both space and time. The vis- 
cous dissipation rate per unit mass is given by 

e (x,t) = K(Vxv) 2 + ^(V-v) 2 ) (1) 

where v is the kinematic viscosity, and we have omitted 
the bulk viscosity, which is negligible for an ideal gas (a 
good approximation for the interstellar medium). The 
divergence term in eq. is necessary for compressible 
flows. Clearly, in the presence of a fluctuating velocity 
field, the dissipation rate is expected to be inhomoge- 
neous. Kolmogorov's 1941 theory uses the average dis- 
sipation rate (as the energy transfer rate in the inertial 
range) and thus implicitly assumes spatial homogeneity 
in the dissipation rate. That theory is therefore not suf- 
ficient to study turbulent heating. Because a significant 
fraction of the kinetic energy is viscously dissipated in 
the finest structures, such as vortex tubes and shocks 



(in supersonic turbulence), which occupy only a small 
volume fraction (or mass fraction for compressible tur- 
bulence), the dissipation rate is strongly intermittent. Its 
probability distribution has a tail much more extended 
than that of a Gaussian distribution, as the result of 
strong local dissipative events. A careful consideration 
of this intermittent distribution is essential for investigat- 
ing turbulent heating. In the following, we briefly review 
the intermittency theory. 

2.1. Intermittency Theory 

Kolmogorov (1962) and Oboukhov (1962) developed 
the first intermittency model for incompressible turbu- 
lence. To account for the fluctuations in the dissipation 
rate, they defined the dissipation rate, e;, at each scale I 
and considered its probability distribution as a function 
of I, which, together with the refined self-similarity hy- 
pothesis (Kolmogorov 1962), was used to predict the scal- 
ing exponents of the velocity structure functions. Gen- 
eralizing to a compressible flow, the dissipation rate per 
unit mass at a scale I is, 

7 TwTTn / p(x + x',t)e(x + x',<)dx' (2) 
pi{x,t)V[l) J 

|x'|<Z 

where V(l) — 4irl 3 /3 is the volume of a spherical region 
of size I and 

MM) = TT^y / P(x + x',t)dx' (3) 

|x'|<Z 

is the average density of that spherical region. From its 
definition, e; is a function of x and we aim to study its 
spatial fluctuations. Note that the average rate, ei, over 
a region of size I defined here is different from vvf /I 2 
used by Falgarone and Puget (1995), which neglects the 
velocity gradient fluctuations in the region, and thus is 
valid only for I smaller than the Kolmogorov dissipation 
scale, T), below which the velocity field is smooth and the 
velocity gradient is approximately constant. 

We point out that the pdf needed to study turbulent 
heating is that at the dissipation scale, the smallest scale 
of fluctuations in the dissipation rate. For this purpose, 
an ideal intermittency model to use is one for the dissi- 
pation range. However, most intermittency models, in- 
cluding the log-Poisson model we will adopt below, are 
for the inertial scales, and we are not aware of theoreti- 
cal models that give the pdf of the dissipation rate in the 
dissipation range. Here, wc take the following approach 
to obtain the pdf at the dissipation scale. We first focus 
on the inertial range and test an intermittency theory for 
the dissipation rate pdf in the inertial scales by compar- 
ing its predicted pdf at each scale with that computed 
from numerical simulations. If the theoretical pdf and 
its scale dependence agree with those from the simula- 
tion data, we extrapolate the model to obtain the pdf 
at the dissipation scale. This assumes that the pdf at 
the dissipation scale can be approximated by that ex- 
trapolated from the inertial scales. Although there is no 
theoretical estimate for its accuracy, the assumption is 
justified by the expectation that the pdf in the dissipa- 
tion range probably connects continuously and smoothly 
to that in the inertial range and is supported by our nu- 
merical results (see below). 
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The probability distribution, P(ei), of the dissipation 
rate, ez, at a scale I can be calculated as 

P(ei) = ^fJ " t))Pi(x, t)dx (4) 

where V is the total volume of the system and p is the 
overall average density in the flow. We have used the 
density as a weighting factor (pi / p, which is unity for in- 
compressible turbulence) to account for the density vari- 
ations in supersonic turbulence. This density or mass 
weighting is supported by Kritsuk et al. (2007), who 
find that replacing the velocity field v by p 1 / 3 v gives a 
3rd order structure function in supersonic turbulence in 
agreement with Kolmogorov's 4/5 law. This is equivalent 
to saying that the density-weighted average dissipation 
rate, (e/), is the same over all the inertial scales. This is 
built in our definition, eq. 1[4"|). 

(ej) = J eiP(ei)dei oc J #(x, <)ej(x,i)dx (5) 

which is independent of I, as follows from eq. @. We 
will therefore denote the average dissipation rate (ej) as 
e. Clearly, e ~ U 3 /L, where U is the rms velocity of the 
flow and L is the integral scale. Another reason for choos- 
ing the density weighted distribution is that we want to 
estimate the mass-weighted average temperature result- 
ing from the turbulent heating of interstellar clouds. 

Kolmogorov (1962) and Oboukhov (1962) proposed 
a log-normal distribution for e;, for which Yaglom 
(1966) gave a justification with a self-similar eddy- 
fragmentation argument. The variance of the log-normal 
distribution is assumed to be proportional to the number 
of cascade steps, oc ln(L/Z), from the integral scale, L, 
to the scale of interest, I. Experiments by Ansemet et 
al. (1984) find that the scaling exponents of the struc- 
ture functions agree with the log-normal model at low 
orders, but depart from it at orders j<b 10. The log- 
normal model has also been shown to violate some theo- 
retical requirements on the structure function exponents 
in incompressible turbulence (Frisch 1995). Numerical 
simulations find that the log-normal distribution of the 
dissipation rate is quite a good approximation in the in- 
ertial range and is valid up to about 4-5 a, beyond which 
a departure from the simulations is clearly seen (e.g., 
Yueng et al. 2006). This departure explains the devia- 
tion of the model from the measured structure functions 
in the experiments. In this work we do not consider 
the log-normal model (which was proposed for incom- 
pressible turbulence), because it gives a poor fit to P{e{) 
computed in our simulations of supersonic turbulence. 

A different intermittency model has recently been pro- 
posed by She and Leveque (1994). The scaling exponents 
of structure functions of order up to 10 predicted by this 
model for incompressible hydrodynamic turbulence agree 
with experimental data with an accuracy of 1%. In their 
original paper, She and Leveque start from the hierarchy 
of dissipation rates at each scale. By invoking a "hidden 
symmetry" that relates the dissipation rate at different 
intensity levels to the strongest dissipative structures at 
a given scale I, they determine the whole hierarchy of 
dissipation rates from the fractal dimension of the most 
intermittent structures and from the scaling of the dis- 
sipation rate of these structures with I. The "hidden 
symmetry" was speculated to be an unknown symmetry 



of the Navier-Stokes equation and was immediately rec- 
ognized as corresponding to a log-Poisson process in a 
multiplicative energy cascade model. 

The model has been successfully applied to incompress- 
ible MHD turbulence (Muller and Biskamp 2000), su- 
personic turbulence (Boldyrev et al. 2002) and super- 
Alvenic turbulent flows (Padoan et al. 2004). These 
studies compute the scaling exponents of the structure 
functions from numerical simulations and compare them 
to the prediction of the log-Poisson model with different 
fractal dimensions for the most intermittent structures. 
They find that a 2D geometry of the most intermittent 
structures, corresponding to current sheets in MHD tur- 
bulence or shocks in highly supersonic turbulence, gives 
an excellent fit to the simulation results. Here, instead of 
considering the scaling exponents of the structure func- 
tions, we focus on the log-Poisson version of the She- 
Leveque model and compare the model directly with the 
probability distribution, P(ei), measured in simulations 
of supersonic MHD turbulence. Our purpose is to check 
whether this model agrees with numerical simulations. 
If so, we can first "calibrate" it by fitting resolved iner- 
tial scales in the simulations, and then extrapolate it to 
high Reynolds number flows. In the following, we give 
a brief description of the log-Poisson model. The inter- 
ested reader is referred to Dubrulle (1994) and to She 
and Waymire (1995) for more details (see also Pan et al. 
2008, which includes both the original presentation by 
She and Leveque and the log-Poisson version) . 

2.2. The Log-Poisson Model 

The log-Poisson model is a multiplicative cascade 
model where the dissipation rates at two scales, l\ and 
I2, are related by a multiplicative factor, Wi t i 2 , 

ei 2 = W hh e h (6) 

Since the average dissipation rate at each scale is equal 
to the overall average dissipation rate e (see above), the 
average of the multiplicative factor, (Wi^), is required 
to be unity. According to the speculation by She and 
Waymire (1995), two types of events determine Wi t i 2 . 
One is the amplification of the dissipation rate that tends 
to cause singular structures at small scales. The other, 
called the modulation-defects by She and Waymire, cor- 
responds to the failure of structures at scale l\ to turn 
into the most dissipative structures at scale I2 in the 
cascade. The first type of events is assumed to give an 
amplification factor of (h/h) 1 , which approaches infin- 
ity as h — > for positive 7. The parameter 7 measures 
how intense the most intermittent structures at scale I2 
can be. Assuming that the kinetic energy per unit mass 
available for dissipation in the most intermittent struc- 
tures is ~ U 2 , where U is the rms velocity at the integral 
scale, and that the dynamical timescale ti in these struc- 
tures follows the regular Kolmogorov scaling, ti oc I 2 / 3 , 
She and Leveque (1994) argued that the energy dissipa- 
tion rate in these intense structures is ~ U 2 /U oc i -2 / 3 , 
which suggests that 7 = 2/3. The modulation-defects 
are assumed to be a discrete process, where each defect 
reduces Wi t i 2 by a factor of (3 and the number of defects, 
n, occurring in the cascade from l\ to I2 is assumed to 
obey a Poisson distribution. The physical meaning of the 
parameter f3 will be discussed below. These arguments 
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Fig. 1. — Logarithm of the dissipation rate (log(e(x)/e), upper 
panel) and density (log(p(x)/p), lower panel) relative to their aver- 
ages on a slice of a snapshot of the simulation. Both the dissipation 
rate and the density have a filamentary structure. Dense filaments 
are usually sites of strong dissipation, but filaments of high dissi- 
pation rate are also found at the interface of low and high density 
regions. As a result, the dissipation rate and the gas density are 
practically uncorrelated (see Fig. l3t, 

are represented by two equations, 

WW = {hlhVP n (7) 

and 

A n 

P(n) = exp(-X hl2 )^ (8) 
n! 

where A;^,, is the average number of defects in the cas- 
cade from I2 to l\. Requiring (Wi t i 2 ) — 1 gives \i x i 2 = 
7^n(7i/Z 2 )/(l — (3). As expected, \i x i 2 is proportional to 
the number of steps in the cascade, ln{l\/l2). 

Taking the logarithm of eqs. ([6]) and (O, we see that the 
distribution of ei is log-Poisson. Given the probability 
distribution at one scale, the distribution of ei at any 
other scale I can be derived using eqs. ©, ([7]) and ([5]). 
The convenient scale to start with is the integral scale, L, 
where the distribution, P^, of the dissipation rate mainly 
depends on the forcing of the flow and thus is probably 
flow-dependent. However, the dissipation rate, e^, at L 
is close to the average, e, and hence Pl is expected to 
be narrow and approximately a delta function. Given 
Pl(Iti€l), we find that 

P(ei)dei = (9) 




Fig. 2. — Logarithm of the projected dissipation rate (upper 
panel) and of the projected density (lower panel) relative to their 
averages, from a snapshot of the simulation. As a result of the lack 
of correlation between the two quantities (see Fig. [3), individual 
structures in the projected dissipation rate do not have an obvious 
counterpart in the projected density. 

00 \ n 

V exp(-\)-rP L (lnfa/e) - jln(L/l) - nln(J3))dln(ei/e) 
* — ' n\ 

where A = \ Ll = jln(L/l)/(l - (3). Note that if P L is 
a delta function, then P(e/) is discrete with a series of 
spikes. 

The physical meaning of the parameter (3 in the log- 
Poisson version of the model is not clear. Deriving the 
moments of the distribution P(ei) from eq. © and com- 
paring with the original version of the model, one finds 
that (3 here corresponds to the same (3 parameter intro- 
duced in the "hidden symmetry" by She and Leveque 
(1994), which has been related to the fractal dimension 
of the most intermittent structures, d, and the parameter 
7 by 7/(1 — (3) = D — d, with D = 3 being the dimen- 
sion of the system. Here we omit the derivation of this 
relation that can be found, e.g., in Pan et al. (2008). 
For vortex tubes, the dimension is d — 1 and for shocks 
d = 2. Because in supersonic turbulence about 1/3 of 
the energy is dissipated in dilatational modes and 2/3 in 
solenoidal modes (Kritsuk et al. 2007), d is expected to 
be between 1 and 2. 

2.3. The Log-Poisson Model for Supersonic Turbulence 
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We compare the log-Poisson distribution, eq. ([9]), with 
results from a numerical simulation of supersonic MHD 
turbulence. We examine how well this model can be 
extended to supersonic MHD flows and what fractal di- 
mension of the most intense dissipative structures gives 
the best fit to the simulation. We take 6 snapshots from 
the 1000 3 Stagger-code MHD simulation in Padoan et 
al. (2007). The simulation adopts periodic boundary 
conditions, isothermal equation of state, random forcing 
in Fourier space at wavenumbers l<fc<2(fc = l cor- 
responds to the computational box size), uniform initial 
density and magnetic field, and random initial velocity 
field with power only at wavenumbers 1 < k < 2. The 
simulation is both supersonic and super-Alfvenic, with 
an rms sonic Mach number of M. s 9, and an initial 
rms Alfvenic Mach number of .M a ,i = 29.7. After a 
few dynamical times, the magnetic energy is amplified 
by the turbulence, and the rms Alfvenic Mach number 
with respect to the rms magnetic field is M a « 2.8, still 
super-Alfvenic. 

In each snapshot, we evaluate the dissipation rate at 
the grid points according to eq. ([1]). A complication 
arises as a result of the limited numerical resolution. The 
dissipation rate should be calculated at the Kolmogorov 
scale using the kinematic viscosity, v ~ c/(na), where 
c is the sound speed, n the number density and a the 
collision cross section. This is not possible because the 
resolution scale in the simulation is much larger than the 
Kolmogorov scale. We cannot use the kinematic viscosity 
to calculate the dissipation rate at the grid points, be- 
cause that neglects the fluctuations in velocity gradients 
below the grid scale and would underestimate the dissi- 
pation rate in each computational cell. To avoid this, we 
are forced to use an effective viscosity at the grid size to 
account for the sub-grid fluctuations erased by the nu- 
merical viscosity. In principle, the effective viscosity at 
each grid depends on the sub-grid fluctuations and (per- 
haps weakly) on the local Kolmogorov scale (which varies 
in space because the kinematic viscosity v depends on 
density) . Due to the lack of information about sub-grid 
scales, we simply assume the effective viscosity is con- 
stant, and normalize it so that the mean dissipation rate 
averaged over all the grid points is equal to the overall 
dissipation rate of the flow, e ~ U 3 /L. With this as- 
sumption, we calculate the dissipation rate in each com- 
putational cell. 

Figure [1] shows the logarithm of the dissipation rate 
(upper panel) and density (lower panel) on a slice of a 
snapshot of the simulation. Both the dissipation rate 
and the density have a filamentary structure. Dense fil- 
aments are usually sites of strong dissipation, but fila- 
ments of high dissipation rate are also found at the in- 
terface of low and high density regions. As a result, the 
dissipation rate and the gas density are practically un- 
corrected. This lack of statistical correlation between 
dissipation rate and gas density is apparent when com- 
paring images of projected dissipation rate and projected 
density, shown in Figure [21 Structures in these two im- 
ages appear to be rather independent of each other. The 
correlation coefficient between dissipation rate and gas 
density, computed over all the computational cells of one 
snapshot, is very low, -0.11, as illustrated by the scatter 
plot of Figure [3l for 10,000 randomly selected computa- 
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Fig. 3. — Dissipation rate versus gas density in 10,000 randomly 
sampled computational cells from one snapshot of the simulation. 
The dissipation rate and the gas density are practically uncorre- 
cted. The correlation coefficient computed over all the computa- 
tional cells of one snapshot is -0.11. 

tional cells of one snapshot. 

From the dissipation rate at each computational cell, 
we compute the average dissipation rate, e;, at any larger 
scale I, using eq. @ (instead of calculating the average 
over spheres, we compute that over cubes of size I from 
the simulation data). We then calculate the density- 
weighted probability distribution, P(ci), of the average 
dissipation rate, e/, over regions of size I, using eq. Q. 
We carry out the same calculation for each snapshot and 
average the distributions from the 6 snapshots. The re- 
sults are plotted in Figure 4 for I = L/64, L/128, L/256, 
and L/512, where L is the size of the simulation box. 
The distribution becomes broader toward smaller scales, 
indicating the dissipation rate is more intermittent at 
smaller scales. Note that the distribution of log(e) is 
strongly skewed toward small values, suggesting that a 
log-normal model cannot give a satisfactory fit to the 
numerical results. 

Figure 4 also shows the probability distribution from 
the log-Poisson model, eq. ([9]), with 7 = 2/3 and (3 = 0.51 
(which corresponds to a fractal dimension of d — 1.64 
for the most intense dissipative structures) at the same 
scales I. As mentioned above, at the box size L, €l = e 
and Pl is close to a delta function. We have set Pl to be 
Gaussian with a small variance that is just large enough 
to make the distribution given by eq. © a smooth func- 
tion. The log-Poisson model with a fractal dimension 
of the most dissipative structures d = 1.64 gives an ex- 
cellent agreement with numerical simulations over more 
than 5 orders of magnitudes in the dissipation rate. The 
agreement for I = L/128 and L/256 is especially remark- 
able, while slight departures for I = L/64 and L/512 are 
perhaps due to, respectively, the memory of large-scale 
motions and the suppression of intermittency by numer- 
ical dissipation. Note that the motivation of the log- 
Poisson model was to explain the tail at large dissipation 
rates corresponding to high-order structure functions. It 
turns out that this model provides a surprisingly good 
fit to the tail of low dissipation rates as well. Padoan et 
al. (2004) obtained a fractal dimension of about 2 from 
MHD simulations with similar Mach numbers as here. In 
that paper, the structure functions were computed with- 
out density weighting, which may explain the difference 
in the resulting dimension of the dissipative structures. 
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Fig. 4. — Density-weighted probability distribution of the dissi- 
pation rate at different scales. The symbols are from our simulation 
of supersonic MHD turbulence, and the curves from the log-Poisson 
model (eq. {9}) with 7 = 2/3 and (3 = 0.51, corresponding to a frac- 
tal dimension d = 1.64 of the most intermittent structures. The 
model gives an excellent fit to the distribution computed from the 
simulation. From bottom to top, the different symbols and the 
curves correspond to the distribution at I//64, L/128, L/256, and 
L/512, where L is the computational box size. The distribution 
becomes broader toward smaller scales, indicating stronger inter- 
mittency at smaller scales. 

Kritsuk et al. (2007) found that two dimensions, d = 1.5 
and d = 2.25, can give acceptable fit to the scaling expo- 
nents for the structure functions of p 1 / 3 v in supersonic 
hydrodynamic turbulence. From their Fig. 2, the fitting 
quality is better for d — 1.5, which is quite close to our 
result here. However, we point out that, except for the 
third order structure function (corresponding to the first 
order moment of the dissipation rate here), the density- 
weighting in J o 1 / 3 v is different from that in eq. (|4|). 

To evaluate the uncertainty in the measured dimension 
d from the average pdf over 6 snapshots, we consider vari- 
ations of the pdf from snapshot to snapshot, which sup- 
posedly include both the intrinsic temporal fluctuations 
of the pdf (if they exist) and the measurement uncer- 
tainty in each snapshot. We think the measurement un- 
certainty at small scales is probably small because of the 
large number of sampling cubes (thus we put more weight 
to the pdf towards smaller scales) and the snapshot-to- 
snapshot variations at these scales may primarily reflect 
the intrinsic temporal fluctuations (which corresponds to 
the variations in the predicted temperature from time to 
time). We calculated the standard deviation of log (P) in 
different e bins over the 6 snapshots. In Fig. 5, we plot 
the 1-cr error bars around the data points for the pdf at 
I = 1/256L, which are normalized by the log-Poisson pdf 
with d — 1.64. Error bars of similar size are also obtained 
for other scales (the size of error bars for the scale, L/64, 
is about 60% larger than those in Fig. 5. This is be- 
cause the measurement uncertainty is larger due to the 
smaller number of sampling cubes.). We see that the 
snapshot-to-snapshot variation is quite small except at 
the far right tail (the effect of which will be discussed 
later in §4.3). From the two curves representing the log- 
Poisson model with d = 1.56 and d = 1.71 (normalized 
to that with d = 1.64), we find that the model gives ac- 
ceptable fits to the measured pdf within the 1-a error 
bars, for a range of values of the dimension d between 
d = 1.56 and d — 1.71. We will evaluate the effect of this 
range of d in our applications. 

We have shown that the log-Poisson model agrees very 
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Fig. 5. — Snapshot-to-snapshot variations of the dissipation rate 
pdf at I = L/256 in the simulation. The data points are the average 
pdf at L/256 (the same asterisk data points in Fig. 4) normalized 
to the log-Poisson pdf eq. J9jl with the best-fit dimension, d = 1.64. 
The error bars indicate the standard deviation (1-cr) of the pdf vari- 
ations over the 6 snapshots. The error bars include both temporal 
variations of the pdf (if they exist) and measurement uncertainties. 
Also shown are the log-Poisson pdfs with d = 1.56 and d = 1.71, 
also normalized to that with d = 1.64. The predicted distribution 
from the model with d in the range (1.56 - 1.71) generally falls 
within the 1 — cr error bars. 

well with the pdf from resolved scales in our simulations 
for supersonic MHD turbulence. As mentioned earlier, 
we are interested the dissipation rate pdf at the Kol- 
mogorov scale, ?y, for applications to turbulent heating, 
and we assumed that the pdf at the scale rj can be ob- 
tained by extrapolating the log-Poisson model, proposed 
for the inertial range only, to the dissipation scale. We 
find that the assumption is well supported by our sim- 
ulation result that pdf in the dissipation scales, such as 
1/512 L, and the resolution scale, is very close to that 
extrapolated from larger, inertial range scales, such as 
1/128 L and 1/256 L. 

Reynolds numbers in interstellar turbulence are much 
larger than can be resolved by current numerical simula- 
tions. Because the intermittency theory can only be ver- 
ified by simulations with limited resolution, the extrap- 
olation of the theory to a realistic value of Re remains 
uncertain. Numerical experiments with much larger reso- 
lution are required to test if the parameters of the theory, 
such as the dimension d, are independent of Re. How- 
ever, based on the usual assumption of turbulence theory 
that the statistics in the inertial range are universal, we 
expect that a theoretical model validated by our simu- 
lation, which resolves at least a small range of inertial 
scales, can be applied to turbulent flows with higher Re. 
If so, the statistics of the dissipation rate extracted from 
the inertial range of our simulation, and the derived di- 
mension d of the log-Poisson model, are universal. 

We will apply the log-Poisson model to different inter- 
stellar clouds using the appropriate values of Re. For 
each type of clouds, we calculate the dissipation scale, 77, 
from the Reynolds number, and use the pdf, eq. Q, at 
I = rj, for the distribution of the turbulent heating rate. 

3. ENERGY BALANCE 

The turbulent heating rate per unit volume is given by, 

T t = n^m/fe(x, t) (10) 

where n and mn are the number density and the mass 
of the hydrogen atom, /i is the mean molecular weight 
(we adopt fi = 2.35 for molecular clouds and /i = 1.23 
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for atomic clouds), and e(x, t) is the turbulent dissipa- 
tion rate per unit mass discussed in § 2. The heating 
rate is spatially inhomogeneous due to the intermittcncy 
of the turbulent dissipation. This gives rise to tempera- 
ture fluctuations. We will use the distribution of e given 
by eq. ([9]), and calculate the mass- weighted average tem- 
perature and the temperature probability distribution in 
dense molecular clouds and in cold diffuse atomic clouds. 

We first consider conventional heating processes in 
these clouds, and then compare their effect with that 
of turbulent heating. Cosmic rays have been thought to 
be the primary heating source for molecular clouds. The 
heating rate by cosmic rays depends on the ionization 
rate, which is very uncertain (e.g., Caselli et al. 1998). 
The measured ionization rate in massive protostellar en- 
velopes, where photoionization is insignificant, is about 
2.6 x 10" 17 s" 1 (Van der Tak and Van Dishoeck 2000). 
Assuming that the energy input per ionization is 20 ev 
(Goldsmith 200f), we take the cosmic-ray heating rate 
per unit volume to be, 

r cr = 0.8 x 10~ 27 n ergs cm" 3 s" 1 (f f ) 

where the number density n is in units of cm -3 . 

In the neutral atomic medium, the dominant heating 
source is the grain photoelectric heating. We use the 
photoelectric heating rate given in Wolfire et al. (2003), 

r pc = 1.3 x 10" 24 ne pc G ergs cm" 3 s" 1 (12) 

where Go is the FUV flux in units of the integrated field 
of Habing (1968) and e pe is the heating efficiency. Fol- 
lowing Wolfire et al., we take Go = 1.7 for the local 
ISM. For the typical cold neutral cloud temperature of 
100 K, density of 50 cm -3 , and electron fraction of 10 -4 
(see Fig. 10 or eq. (C15) in Wolfire et al. 2003), we get 
e pc = 1.4 x 10" 2 (using eq. (20) in Wolfire et al. 2003), 
and hence r pc = 3.1 x 10" 26 n ergs cm -3 s" 1 . 

For molecular clouds, we take the cooling rate per unit 
volume, A g , from Table 2 of Goldsmith (2001), which 
gives coefficients and indices for power-law fits to the 
temperature dependence of the cooling rates at densi- 
ties from 10 2 to 10 7 cm -3 . We use a linear interpola- 
tion to get coefficients and indices for densities not given 
in the table. The table is only for velocity gradient of 
1 km s" 1 pc _1 . According to eqs (2), (3) and (4) in 
Goldsmith and Langer (1976), the cooling rate per par- 
ticle depends on the velocity gradient, dv / dr, and on the 
density only through their product, n/(dv/dr). There- 
fore, Table 2 of Goldsmith (2001) can be converted to one 
for the cooling rate per particle as a function of n{dv/dr). 
We then obtain the cooling rate per particle for a given 
density and velocity gradient in the range of n/ (dv / dr) 
corresponding to the density range in the table. 

For cold neutral HI clouds, the dominant cooling pro- 
cess is the line cooling by the CII 158 /im fine structure 
transition. Considering the low electron fraction (10" 4 ) 
in these clouds, we neglect the contribution to the CII 
cooling from collisional excitation by electrons. In § 4.3, 
we compute the fraction of gas in cold diffuse clouds 
with high enough temperature ( ^ 1000 K) to activate 
the chemical reactions producing CH + . We therefore in- 
clude cooling by the OI 63 fim line, which becomes more 
important than CII cooling at T > 1000 K. We adopt 
the cooling rates for CII and OI lines from Wolfire et al 
(2003). Note that this cooling rate is smaller than that 



given by Dalgarno and McCray (1972) by a factor of a 
few. This is because Wolfire et al. (2003) adopted C and 
O abundances derived from observations of UV absorp- 
tion lines by the interstellar gas over many lines of sight 
studied by HST/FUSE, which are smaller than the solar 
abundances used in Dalgarno and McCray (1972) by a 
factor of 2-3. 

In our calculations, we will consider energy exchange 
between gas and dust grains and solve the coupled energy 
balance equations for gas and dust grains simultaneously. 
For the energy transfer rate from gas to dust grains, A gc j, 
we take the formula, 

A gd = 0.7 x 10" 33 n 2 T 1/2 (T - T d ) ergs cm" 3 s" 1 (13) 

from Goldsmith (2001), where T and Td are, respectively, 
the gas and the grain temperatures in units of K (see also 
Black 1987, where the rate is larger by a factor of a few). 
Eq. (fl3|) assumes a dust to gas ratio of 0.01, grain size 
of 1.7 x 10" 5 cm, grain density of 2 g cm" 3 , and an 
accommodation coefficient of 0.3 for grains and H2. If 
T > Td, the gas heats the dust and A g d is a cooling term 
for the gas; vice versa, if T < Td, A g d is a heating term 
for the gas. 

Dust grains in molecular clouds are heated by the 
diffuse UV-visible-IR interstellar radiation. Following 
Goldsmith (2001), we adopt a flux of 5.3 x 10 -3 ergs cm 
" 2 s" 1 for this radiation field. Using the above values for 
the dust to gas ratio and the grain cross section (calcu- 
lated from the grain size and density), the dust heating 
rate by the diffuse interstellar radiation is, 

r d = 3.9 x 10" 24 nx ergs cm" 3 s" 1 (14) 

To account for the attenuation of the radiation field due 
to dust extinction, we have reduced the heating rate by 
a factor of \. We take x to be ~ 10" 1 — 10" 2 , corre- 
sponding to a visual distinction of 2-5 mag., typical of 
a molecular cloud of size 1 pc. The value of x affects 
only the dust temperature in such a cloud (§4.1). For in- 
ner regions of high-mass star-forming cores, the radiation 
from newly formed stars is the primary heating source of 
dust grains (see § 4.2). The heating rate from the stellar 
radiation gives a high dust temperature. In that case, 
collisions of gas molecules with dust grains could provide 
a significant energy source for the molecular gas through 
A g d- For dust heating in atomic clouds, we use eq. (|14| 
with x = 1- 

The dust cooling rate is calculated by integrating over 
the wavelength, A, the absorption efficiency, Q\, times 
the Planck function (Hollcnbach and McKee 1979). For 
A < 100/im (which corresponds to Td 140 K), the 
absorption efficiency decays as A" 2 (Draine and Lee 1984; 
Ossenkolf and Henning 1994). Using Q x ^ 4.5 x 10" 5 at 
the reference wavelength of A = 790/im (Goldsmith 2001; 
Ossenkolf and Henning 1994), a numerical evaluation of 
the integral gives, 

A d = 4 x 10" 31 nT d 6 ergs cm" 3 s" 1 (15) 

(Note that this result of our numerical calculation is 60 
times larger than that given in Goldsmith (2001)). In 
our calculations below, the dust temperature is smaller 
than 140 K, eq. (fl~5| is thus valid for the temperature 
range of interest here. 

Putting together heating and cooling for gas and dust 
grains and their energy exchange, the energy balance 



equations for gas and dust grains are, respectively, 

r t + r cr + r pc - A g - A gd = o (16) 

and 

r d + A gd - A d = (17) 

The energy balance assumption is here justified because 
the cooling timescale of interstellar clouds we consider in 
the next section (~ 10 4 yr) is smaller than their dynami- 
cal timescale (~ 10 6 yr). We will solve these two coupled 
equations numerically. 

4. APPLICATIONS 

We apply the intermittent turbulent heating model to 
three different interstellar medium environments: Molec- 
ular clouds, high-mass star-forming cores, and cold dif- 
fuse HI clouds. The spatial variations in the turbulent 
heating rate, I\, given by eq. (fT0|) . may result in large 
spatial fluctuations of the gas temperature. Given the 
distribution function of T t , we can calculate the temper- 
ature probability function by solving eqs. (fl"6|) and (JTTJ) . 
However, there is a complication in converting the mass- 
weighted distribution of the heating rate to the mass- 
weighted temperature distribution. The various terms 
in eqs. (fl6| and (|T7|) have different density dependence 
and thus an exact calculation of the mass-weighted tem- 
perature distribution requires the joint statistics of den- 
sity and dissipation rate fluctuations, which is beyond 
the scope of this study. We will therefore solve the en- 
ergy balance equations using simply the average gas den- 
sity. This is equivalent to assuming that the density field 
and the dissipation field are not correlated, which is sup- 
ported by their low correlation coefficient of -0.11 and 
by Figures [2] and [3l Although the temperature distribu- 
tion is not mass-weighted in an exact way, the resulting 
temperature distribution and the average temperature 
from our calculations are mass-weighted in the sense that 
they are computed from the mass-weighted distribution 
of the dissipation rate. In the case of molecular clouds, 
the mass-weighting of the temperature distribution dis- 
cussed in § 4.1 is almost exact. The cooling rate in these 
clouds has an almost linear dependence on density (Gold- 
smith and Langer 1978) so, except for the negligible cou- 
pling term A gd , the density dependence in all the terms 
of eq. (fl6|) is the same. 

We point out that the nearly independence between 
the dissipation rate (per unit mass or per particle) and 
the density does not imply that the temperature and the 
density are decoupled. The gas temperature is calculated 
from the balance between heating and cooling. The cool- 
ing rate per particle depends on the density and thus 
the resulting temperature is correlated with the density 
through the cooling rate. In this sense, the correlation is 
similar to the case with a uniform heating rate (as usually 
assumed for traditional heating sources, e.g., the cosmic- 
ray heating): denser regions tend to have a lower temper- 
ature than lower density ones. However, the intermittent 
fluctuations in the turbulent heating rate may also give 
rise to significant scatter around this general trend. This 
scatter may be probed through high-resolution molecular 
emission line maps. In such maps, the observed scatter 
around the predicted correlation between the projected 
temperature and density (i.e., averaged along the line of 
sight) would be smaller than that in 3D space, as can be 
seen by comparing Fig. 1 with Fig. 2. 



4.1. Molecular Clouds 

We assume the following Larson's relations for the ID 
velocity dispersion and the number density in molecular 
clouds, 

cr v = l(L/l pc) 04 km s~ 4 (18) 

and 

n = 2000(L/1 pc)" 1 cm" 3 (19) 

where L is the size of the cloud. The average dis- 
sipation rate per unit mass is e = l/2(v / 3cr v ) 3 /_L = 
0.84 x 10~ 3 (L/1 pc) 2 ergs g -1 s _1 , where we used the 
conversion from the ID rms velocity, <r v , to the 3D rms 
velocity, U, i.e., U — \/3a v . We have assumed that the 
turbulence is driven at the length scale of the cloud (see 
Basu and Murali 2001, who find that driving at smaller 
scale, for example by stellar outflows, would result in 
a turbulent heating rate that gives a CO luminosity in 
excess of the observations). Therefore, the average tur- 
bulent heating rate per unit volume, given by eq. (|10[) . 
is, 

f t = 3 x lQ- 27 n(L/l pc) 2 ergs cm" 3 s" 1 (20) 

In reality, the turbulence may be driven at a scale much 
larger than the cloud size, as suggested by various stud- 
ies of the velocity scaling in molecular clouds (e.g., Os- 
senkopf and Mac Low 2002; Heyer and Brunt 2004; 
Padoan et al. 2006). However, from eq. (|20p. T t is al- 
most independent of L, and therefore eq. (|2t)]) is accurate 
within a factor of a few. By adopting the cloud size as 
the driving scale, we make a conservative assumption, 
which may slightly underestimate the turbulent heating 
rate and the resulting gas temperature. 

Comparing eq. (|20|) with eq. (fTTj) . we see that the av- 
erage turbulent heating rate in molecular clouds is of the 
same order of the cosmic ray heating rate, and possi- 
bly a few times larger. This average turbulent heating 
alone could maintain a temperature of ~ 10 K in molec- 
ular clouds, and thus the cosmic ray heating, whose rate 
is highly uncertain, is not even needed to explain the 
observed gas temperatures. Neglecting the spatial fluc- 
tuations in the turbulent heating rate, and using the av- 
erage heating rate from eq. (|20|) in eq. (fl"6]) , and the dust 
heating rate from eq. (TT4|) in eq. (jTTJ) , we find that the 
turbulent heating alone gives a temperature of 17 K for a 
1-pc cloud. Similar temperatures are obtained for clouds 
of size 0.1 pc and 10 pc, because the heating rate from 
eq. ([20]) depends very weakly on I. In the density range 
of these clouds (200-2 xlO 4 ), the coupling between gas 
and dust is week, so A gd is negligible in comparison with 
the gas and dust heating rates, and can be neglected. 
The dust and the gas temperatures can be obtained sim- 
ply from their own energy balance (and thus the gas tem- 
perature is independent of x m ec l- CZIOj an d, from the 
interstellar UV-visible-IR radiation field, the dust grains 
achieve a temperature of 7 K if \ — 10~ 2 , or 10 K if 

x = io- 1 . 

A realistic calculation for the gas temperature needs to 
take into account the intermittent distribution of the tur- 
bulent heating rate. Note that the cooling rate in molec- 
ular clouds depends on temperature quite sensitively; for 
example at n = 1000 cm~ 3 , A g oc T 2A (Goldsmith 2001). 
Therefore, given the average temperature, the presence of 
temperature fluctuations makes the overall cooling rate 
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larger. Conversely, given the average heating rate (and 
thus the overall cooling rate), the presence of heating 
rate fluctuations (which give rise to temperature fluctu- 
ations) leads to an average temperature lower than from 
uniform heating. The larger the fluctuation amplitude 
(for example from stronger intermittency in the heating 
rate), the smaller the average temperature. This effect 
can be significant because the log-Poisson distribution of 
the heating rate, eq. ©, is very skewed toward low values 
(see Figure^). In other words, the turbulent heating rate 
is below its average value in a large fraction of the mass. 
Because the hotter gas cools faster than the cooler gas, 
the contribution to the average temperature is mainly 
from the low-dissipation-rate part of the distribution. 

To quantify this intermittency effect, we assume the 
distribution of the turbulent heating rate, T t , follows 
the log-Poisson distribution, eq. ©, which we evalu- 
ate at the dissipation scale, rj. The dissipation scale is 
77 = Li?e -3 / 4 , where the Reynolds number, Re, is defined 
as Re = UL/v. The kinematic viscosity is calculated as 
v = c/(na) (see § 2). We adopt a typical value for the 
cross section, a ~ 10~ 15 cm 2 , and use the average tem- 
perature to calculate the sound speed, c. By numerically 
solving eq. (fT6|) and eq. (fT7|) . we convert the distribution 
of T t into the temperature distribution. We assume the 
dimension of the intermittent dissipative structures to be 
d = 1.64. The cumulative probability of the gas temper- 
ature is plotted in Figure [6] for a 1-pc cloud (solid line). 
We find that the average temperature of the distribution 
is 8.5 K. As expected, this is lower than 17 K, obtained 
using the average dissipation rate, but still a reasonable 
mean temperature for molecular clouds, in the absence 
of cosmic ray heating. 

We estimate the effect of the uncertainty in d measured 
from snapshot-to-snapshot variations, 1.56 < d < 1.71, 
and find that this range of values of d corresponds to 
a range of 8.0-8.9 K for the average temperature, i.e., 
an uncertainty of only ~ 6% around the value of 8.5 K, 
corresponding to d = 1.64. The effect on the average 
temperature is so small because the average tempera- 
ture is ~ (e 1 / 2-1 / 3 ) (since the cooling rate is ~ T 2 ~ 3 for 
the density, 10 3 — 10 4 cm" 3 , in the molecular clouds of 
interest here), which is mainly contributed by the part 
of the pdf between its peak and the average dissipation 
rate. From the error bars in Fig. 5, one can see there is 
little snapshot-to-snapshot variation in this part of the 
pdf. To further illustrate the effect of the degree of in- 
termittency, we also calculated the average temperature 
using different values of the dimension of the dissipa- 
tive structures. For example, if d = 1, as in the case 
of incompressible flows, the dissipation rate distribution 
would be less intermittent and the average temperature 
would be 11 K. On the other hand if d = 2, the average 
temperature would be lower, about 6 K. 

Although the turbulent heating alone gives a reason- 
able average molecular cloud temperature of ~ 8.5 K, 
a significant mass fraction of the clouds receives little 
heating from turbulent dissipation, because the distri- 
bution of the dissipation rate is skewed toward small 
values. For d = 1.64, about 20% of the mass is not 
significantly heated by turbulence and, in the absence 
of other heating processes, this mass fraction would be 
at the same temperature of 2.7 K as the cosmic mi- 




1 10 100 

T(K) 

Fig. 6. — Cumulative probability of temperature smaller than T 
in a molecular cloud of size 1 pc. The three curves correspond to 
the results for different cosmic ray heating rates. Without cosmic 
ray heating (solid line), the average temperature from turbulent 
heating is 8.5 K, and 20% of the gas mass is essentially not heated, 
due to the intermittency of the turbulent heating rate. When in- 
cluded (dashed and dotted lines), the cosmic ray heating sets a 
lower limit for the temperature of 9.5 K and 14 K (corresponding 
to the uniform temperature in the absence of turbulent heating), 
for the two heating rates considered here. In these two cases, the 
average temperature is 13 K and 17 K. 

crowave background (see Figure [6]). The cosmic ray 
heating would of course dominate in regions with es- 
sentially no turbulent heating. In Figure [TJ] we also 
show the temperature distribution including both tur- 
bulent and cosmic ray heating. To account for the un- 
certainty in the cosmic ray heating rate, we use two 
values, r cr = 0.8 x 10~ 27 n ergs cm -3 s _1 (eq. ([TTjl : 
dashed line) and r cr = 2.0 x 10~ 27 n ergs cm -3 s _1 
(dotted line). The heating rate by cosmic rays sets the 
lower temperature limit (corresponding to the uniform 
temperature in the absence of turbulent heating). For 
T = 0.8 x 10~ 27 ergs cm -3 s _1 , this lower limit is 9.5 K. 
In this case, the average temperature is 13 K and most 
gas (80% of the mass) has a temperature lower than the 
average. The larger cosmic ray heating rate chosen here, 
r cr = 2.0 x 10~ 27 n ergs cm -3 s" 1 , gives a lower limit of 
14 K and an average temperature of 17 K. The tail of 
large turbulent dissipation rate gives a small mass frac- 
tion of gas with high temperature. We find that only 1% 
of the gas could be hotter than 40 K. 

We have shown that the turbulent heating by itself 
can maintain an average temperature of about 8.5 K in 
molecular clouds. Due to the intermittency of this heat- 
ing process, a finite but small fraction of gas has temper- 
ature higher than 40 K, and a significant fraction of the 
gas is essentially not heated by turbulence, but rather 
by cosmic rays. In the presence of cosmic rays, turbu- 
lent heating increases the average temperature by a few 
degrees, from 9.5 K to 13 K for the cosmic ray heating 
rate given by eq. (fTT|) . 

Similar calculations can be done for dense dark cloud 
cores. We use the dense cores in the Orion A molecu- 
lar cloud complex as an example (e.g., Li et al. 2003, 
2007), and focus on the ones with no apparent internal 
sources for simplicity 1 . The sizes of these cores are typ- 
ically 0.05-0.2 pc (e.g., Tatematsu et al. 1993, Li et al. 

In dense cores containing young stellar objects (e.g., those dis- 
cussed in Johnstone and Bally 2006), one needs to consider heating 
from the central sources and the calculation would be similar to 
that for the high-mass star-forming cores to be discussed in §4.2 
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2003, Ikeda et al. 2007) and the typical average density 
is 10 4 — 10 6 cm~ 3 . Observations using tracers such as 
CS (e.g., Tatematsu et al. 1993), NH 3 (e.g., Li et al. 
2003), H 13 CO+ (Ikeda et al. 2007), and other molecu- 
lar transitions show that the spectral lines have FWHM 
(full width at half maximum) in the range 0.5-2 km s , 
corresponding to 0.22-0.89 km s _1 for the ID rms veloc- 
ity, cr v . We emphasize that we evaluate the dissipation 
rate based on the rms velocity, while observational papers 
usually give the FWHM. We obtain the ID rms velocity 
from the conversion cr v = FWHM/2.355. 

We find that the typical values of cr v = 0.5 km s _1 and 
L = 0.1 pc follow quite closely Larson's relation given 
in eq. (|18[) . implying that the turbulent heating rate is 
similar to that in typical regions of molecular clouds dis- 
cussed above. A complication in evaluating the temper- 
atures in these cores is the possible presence of external 
UV sources (e.g., from the Trapezium star cluster, see 
Li et al. 2003). We also need to consider the inner and 
outer regions of these cores separately. The central re- 
gion is self-shielded from the diffuse interstellar radiation 
field or the external UV sources, and thus has a lower 
dust temperature than the outer regions if there are no 
internal sources. Furthermore, the gas and dust temper- 
atures are well coupled in the central region, due to the 
high density (~ 10 6 cm -3 ) (see §4.2 for a more detailed 
discussion on the gas-dust coupling). We find that in- 
cluding turbulent heating has little effect to the coupled 
gas and dust temperature in the central region and the 
temperature there is essentially determined by the dust 
heating rate from the diffuse interstellar radiation field 
and the external UV sources. 

At the lower density outer edge, gas and dust are ther- 
mally decoupled, and we find that the effect of turbulent 
heating to the gas temperature is similar to that in the 
molecular clouds discussed above. As an example, we 
assume that the dust heating is from the diffuse inter- 
stellar radiation field and the rate is given by eq. (TT4"|) . 
with x ~ 1 at the edge. This rate gives a dust tempera- 
ture of 15 K, close to the average dust temperature of 17 
K found in the observations by Lis et al. (1998). Assum- 
ing the gas density is 10 4 cm -3 at the edge, our model 
with d = 1.64 shows that turbulent dissipation alone can 
heat the gas to an average temperature of 11 K, close 
to the temperature of 10 K due to cosmic-ray heating 
alone. Combining the comic-ray heating with the turbu- 
lent heating results in a gas temperature of 14 K, which 
is in general agreement with the observed temperature 
(e.g, Li et al. 2003). The effect of turbulent heating 
here is thus the same as that in molecular clouds stud- 
ied above. Some cases of higher gas temperature have 
also been found in dense cores in Orion, which can be 
explained by external UV sources (or internal sources if 
they are present). UV sources directly heats dust grains 
and increase the dust temperature. Since in the denser 
inner regions gas and dust are thermally coupled, the 
gas temperature there may be increased to the observed 
higher values by UV sources through energy transfer with 
dust grains. (The gas temperature at the outer edge with 
a density ~ 10 -4 cm -3 may not increase considerably 
because the gas and dust temperatures are there decou- 
pled.) 

In summary, turbulent heating in central regions of the 
dense cores in Orion have negligible effects to both gas 



and dust temperatures, while its role for the outer re- 
gions is similar to that in the molecular clouds following 
Larson's relation. Turbulent heating alone can heat the 
gas to a temperature of 11 K, and a few degrees higher in 
combination with cosmic ray heating. The effect of tur- 
bulent heating is more prominent in objects with much 
larger turbulent intensity, such as in the high-mass star- 
forming cores to be discussed in § 4.2, whose average 
temperature is significantly increased, or in places where 
the intermittent tail of the heating rate pdf plays an im- 
portant role for the chemistry, such as in HI clouds, which 
we study in § 4.3. 

4.2. High-Mass Star-Forming Cores 

We now apply the turbulent heating model to the high- 
mass star-forming cores observed in various CS transi- 
tions by Plume et al. (1997). These cores were originally 
selected by the presence of water masers, suggesting the 
formation of massive stars, as also shown by their large 
FIR luminosity (Mueller et al. 2002). The characteristic 
mean density of these cores is 10 6 cm -3 , the mean size 
0.3 pc (average values over the objects listed in Plume 
et al. (1997) and Shirley et al. (2003) ). These cores are 
very massive, with a typical mass of ~ lOOOAf© (Plume 
et al. 1997, Shirley et al. 2003) (much more massive 
than the dense cores in Orion discussed in §4.1). We 
are interested in these cores because of their extremely 
large turbulent intensity. The observed FWHM is in 
the range of 2-12 km s _1 (Shirley et al. 2003) with 
an average of ~ 6 km s^ 1 . This average correspond to 
a ID rms velocity of 2.5 km s _1 . The average turbu- 
lent dissipation rate in these dense cores is very large, 
e = 4.4 x 10~ 2 ergs g _1 s _1 , approximately 50 times 
larger than that in a molecular cloud following Larson's 
relations (see § 4.1). However, this heating rate per unit 
mass, which can be written as ~ 0.02 Lq/Mq, is much 
smaller than the average (FIR) bolometric luminosity to 
mass ratio, 140 L©/M , found by Mueller et al. (2002). 
Clearly, turbulent heating in these cores could not ac- 
count for even a tiny fraction of the FIR radiation. The 
FIR emission must be primarily due to the processed UV 
photons from the central massive stars, which directly 
heat the dust grains, resulting in a fairly high dust tem- 
perature. The primary gas heating source in these cores 
is thought to be the energy exchange with dust grains 
by collisions (Evans 1999). We want to investigate if the 
turbulent heating may provide a significant energy source 
for the gas in these dense cores and to estimate the effect 
of the turbulent heating on the gas temperature. 

The gas density and the dust heating rate from the cen- 
tral stellar sources (which can be estimated from the dust 
temperature) are needed to calculate the gas tempera- 
ture from eqs. (fTE)) and (fT7|) . The gas density, through 
the collision frequency, plays a crucial role in determin- 
ing how well gas and dust grains are thermally coupled 
(see eq. (TT3")) ). The degree of coupling decreases with 
the distance to the center, r, as the density decreases. 
The density profile in these cores is usually assumed to 
be a power-law, n(r) oc r~ Q , where a is in the range 
0.5-2.5 (Mueller et al. 2002). In this range of a, most 
of the gas mass is in the outer regions at large r; even 
for the steepest slope, a — 2.5, about 70% of the gas 
mass is outside r = 0.1 pc, for a core of size 0.3 pc. 
Since a large fraction of the gas mass is in the outer re- 
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gions, one may expect that the density there is about 
the average density, 10 6 cm -3 , given earlier. However, 
the measured average density of 10 6 cm" 3 is evaluated 
from radiative transfer modelling of the emission lines of 
tracer molecules (Shirley et al. 2003), and thus proba- 
bly reflects the average density along the line of sight, 
which is dominated by inner regions for a density profile 
steeper than r . In that case, the density at the outer 
regions is much smaller than 10 6 cm -3 . 2 For example, 
if a = 2, normalizing the line-of-sight average density to 
10 6 cm -3 and assuming that the radius at the inner edge 
of the core is 100 AU, we find that the density at 0.1 pc 
is 10 4 cm" 3 (while the density at 0.01 pc is close to the 
measured average density of 10 6 cm" 3 ). To account for 
the possibility of a flatter density profile, we take n at 
0.1 pc to be 10 4 — 10 5 cm' 3 . We did similar estimates 
for the density at inner regions and found that, almost 
independent of the slope, a, in the range 0.5-2.5, the gas 
density at 0.01 pc is always close to the measured average 
density of 10 6 cm -3 . 

Our calculations show that 10 6 cm" 3 is a critical den- 
sity for the gas and dust temperatures to be well coupled. 
Thus the gas temperature behavior within 0.01 pc is dif- 
ferent from that outside 0.01 pc. In our calculations, 
we will refer to two characteristic values of the radius, 
r = 0.01 pc and r = 0.1 pc, representative of the dif- 
ferent gas temperature behaviors in the inner and outer 
regions of these cores. 

The dust temperature depends on the heating rate 
from the central stellar source, which decreases with in- 
creasing radius roughly as r~ 2 . Using the cooling rate 
given by eq. (fl"5"j). we obtain Td oc r -1 / 3 , which is close 
to the power-law fit (Td oc r" 0,4 ) to the numerical result 
from the radiative transfer model in Mueller et al. (2002) 
(the difference in the power-law indices here may suggest 
that a radiative transfer calculation gives a dust heating 
rate decreasing faster than r -2 , or a dust cooling rate in- 
creasing with Td more slowly than oc T^). We find that 
the energy transfer from dust to gas or from gas to dust 
(which is possible in the presence of turbulent heating of 
the gas) is negligible in comparison with the dust heating 
rate from the stellar photons, so the dust temperature is 
completely determined by the central stellar source. For 
a typical bolometric luminosity of 10 4 T©, the dust tem- 
perature decreases from about 65 K at r = 0.01 pc, to 
about 25 K at r = 0.1 pc (see Fig. 9 in Mueller et al. 
(2002)). 

We now calculate the gas temperatures at 0.01 pc and 
0.1 pc. We first study the inner region using Td = 65 K, 
and a density of 10 6 cm -3 . This dust temperature im- 
plies a heating rate of 3 x 10" 14 ergs cm" 3 s" 1 by UV 
photons from the central stars (estimated from the dust 
cooling rate, eq. (15])). Fixing this dust heating rate, we 
first calculate the gas temperature assuming the gas is 
heated only by collisions with dust grains. By solving 
eqs. (|16|) and (fT"7|) . we find a gas temperature of 57 K, 

2 Plume et al. (1997) found that the average density, ~ 10 6 
cm -3 , derived from their radiative transfer model, gave a core 
mass larger than the virial mass by more than an order of magni- 
tude. They realized that the mass discrepancy might have origi- 
nated from the assumption of uniform density in their model. This 
suggests that the true average density, determined mainly by the 
outer regions, is smaller than 10 6 cm -3 by one or two orders of 
magnitude, supporting our choice of 10 4 -10 5 cm . 



suggesting a strong thermal coupling between gas and 
dust, as mentioned above. The gas heating rate per unit 
volume by dust grains is 3.7x 10 -20 ergs cm" 3 s" 1 , which 
is only ~ 10~ 5 times the heating rate of the dust grains 
by the stellar radiation field. Thus the dust temperature 
is not affected by the energy transfer to the gas. We 
then calculate the gas temperature at r = 0.01 pc as- 
suming the gas is heated both by the dust and by the 
turbulence. The average turbulent heating rate in these 
cores is 1.3 x 10 -19 ergs cm -3 s" 1 , which is larger than 
the rate from collisions with dust grains when the tur- 
bulence is neglected. The inclusion of turbulent heating 
should therefore result in a gas temperature larger than 
57 K. Using again the intermittency model with d = 1.64, 
we obtain a mean gas temperature of 64 K, very close to 
the dust temperature of 65 K. In summary, the gas and 
the dust temperatures in this dense inner region are very 
well coupled, but the turbulent heating can still increase 
the gas temperature by a few degrees. 

The outer region is interesting because it contains most 
of the mass of the core and because its density is lower 
and so gas and dust are less well coupled. We carry out 
the same calculation as above for the outer region, at 
r = 0.1 pc, assuming a density of n = 10 4 -10 5 cm -3 and 
a dust temperature of Td = 25 K. Neglecting turbulent 
heating, we obtain a gas temperature of 17 K and 6 K, 
for n — 10 5 and 10 4 cm -3 respectively (including cos- 
mic ray-heating increases the temperature to 18 K and 
10 K respectively). In this case the gas temperature is 
considerably smaller than the dust temperature due to 
the weak coupling (while with n = 10 6 cm" 3 the gas 
and dust temperatures would be almost equal). As we 
did for the inner region, we now include the turbulent 
heating. For illustration, we assume the average dissi- 
pation rate is independent of radius. Without energy 
transfer between dust and gas, we find that the intermit- 
tent turbulent heating alone would give a temperature of 
36 K and 35 K, for n = 10 4 and 10 5 cm" 3 respectively. 
The gas temperature is here nearly independent of the 
density because at densities above 10 4 cm -3 the cooling 
rate per particle is almost constant (see Goldsmith and 
Langer 1978). Because in this outer region the thermal 
coupling between dust and gas is weak, we find that in- 
cluding both the energy transfer between dust and gas 
and the intermittent turbulent heating, the gas temper- 
ature is still 35-36 K. In other words, in the outer region 
of the cores, containing most of the core mass, the gas is 
thermally decoupled from the dust and its temperature 
is completely determined by the turbulent heating. Sim- 
ilar to the case of molecular clouds, we find that range 
of d from 1.56 to 1.71 results in a small range for the av- 
erage temperature, within 1.5-2 K (i.e., ~ 5%), around 
the value 35-36 K obtained for d — 1.64. 

Our result is in general agreement with the observa- 
tions of high-mass protostellar objects by Leurini et al. 
(2007). By fitting the observed methanol spectral lines, 
they derive the physical parameters of a model consist- 
ing of an inner core, an extended component, and an 
outflow. In particular, they find the kinetic temperature 
is in range of 22-40 K in the extended component (cor- 
responding to the outer regions in our model) where the 
FWHM linewidth is 2.8-4.2 km - 1 (see their Table 5) 
over an average size of 0.1 pc (smaller then the value of 
0.3 pc we adopted above). With these specific parame- 



12 



ters, and the density range of 10 5 — 10 6 cm" 3 (Table 5 in 
Leurini et al. 2007), our model would predict a temper- 
ature range of 28-39 K. Although the velocity dispersion 
here is lower than adopted in our example, the size is also 
smaller, resulting in a dissipation rate and a temperature 
comparable to those in our example above. Future obser- 
vations with high spatial resolution will provide a better 
test of our model. 

We conclude that, although negligible in comparison 
with the dust heating rate by the central stellar sources, 
turbulent heating in high-mass star-forming cores pro- 
vides an important energy source for the gas. For a char- 
acteristic density of 10 6 cm" 3 in the inner regions, the 
gas and dust temperatures are well coupled. The tur- 
bulent heating increases the gas temperature only by a 
few degrees. This slight temperature increase in the in- 
ner regions may not have important observational or dy- 
namical consequences. On the other hand, the turbulent 
heating has significant effects in the outer regions, where 
the density is lower and the coupling of dust and gas is 
weaker than in the inner regions. For a characteristic 
dust temperature of 20-30 K and a density of 10" 4 -10" D , 
the turbulent heating increases the gas temperature from 
less than 10-20 K (lower than the dust temperature) up 
to about 36 K. This large increase in the gas tempera- 
ture by a factor of 2-3 due to turbulent heating, in the 
region containing most of the core mass, may have im- 
portant effects on the evolution of the cores, on their 
star formation efficiency, and on their stellar initial mass 
function. Future observations of these cores may help de- 
termine their gas temperature with sufficient accuracy to 
constrain the relative importance of turbulent heating. 

4.3. Cold Diffuse HI Clouds 

The observed abundance of CH + , HCO + , and OH 
molecules in cold diffuse clouds (Crane et al. 1995, Gre- 
del 1997, Lucas and Liszt 1996, Liszt and Lucas 1996) 
suggests the existence of hot regions inside these clouds. 
The activation temperature of reactions producing these 
molecules is > 1000 K, while the cloud average temper- 
ature is < 100 K. Dissipative heating by MHD shocks 
and vortex tubes have been proposed to explain the exis- 
tence of hot regions (see, e.g., Pety and Falgarone 2000). 
These dissipative structures are the most intermittent 
structures in the log-Poisson model of the turbulent en- 
ergy dissipation. Since the model gives the distribution of 
the dissipation rate at all intensity levels, it can be used 
to calculate the cumulative probability corresponding to 
the mass fraction of regions with T > 1000 K. Falgarone 
and Puget (1995) performed a similar calculation. They 
used the probability distribution of the velocity differ- 
ence from experiments with incompressible flows, where 
the strong dissipative structures responsible for the hot 
regions are vortex tubes. It is not clear whether their 
result applies to supersonic turbulence in diffuse neutral 
clouds, where both shocks and vortex tubes contribute to 
high temperatures. The log-Poisson intermittency model 
for supersonic turbulence includes both types of dissipa- 
tion structures. We will compare the result by Falgarone 
and Puget (1995) with that from the log-Poisson model 
for incompressible turbulence and for supersonic flows 
presented in § 2. 

For the characteristic density and length scale of cold 
neutral clouds we adopt n = 50 cm" 3 and L=l pc (Helies 
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Fig. 7. — Cumulative probability of regions with temperature 
smaller than T in diffuse HI clouds. The three curves correspond 
to the results for three different values of the rms turbulent velocity, 
<7 V = 1, 2, and 3 km s , at the scale of 1 pc. The mass fraction of 
gas with T > 1000 K for these three cases is 3.6X10 -5 , 3 X 10" 3 , 
and 0.02, respectively. The observed CH + abundance requires that 
ov is at least 3 km s" 1 at 1 pc. Most of the gas has a temperature 
below 100 K. 

and Troland 2003). The observed ID rms turbulent ve- 
locity, ov, at the scale of 1 pc, is between 1 and 3 km s" 1 
(Hcnnebelle et al. 2007). The average turbulent heat- 
ing rate is 1.5 x 10~ 27 (er v /l km s -1 ) 3 ergs s" 1 per H 
atom. We will treat a v as a parameter and study the 
dependence of the mass fraction of hot regions on a v . As 
discussed in § 3, the photoelectric heating rate in cold 
neutral clouds is 3.1 x 10 -26 ergs s _1 per H atom. Ne- 
glecting turbulent heating, this gives a temperature of 53 
K using the CII and OI cooling rates from Wolfire et al. 
(2003). At n = 50 cm" 3 , the gas-dust coupling is very 
weak and their energy balance equations can be solved 
separately. 

We now consider the effect of turbulent heating. First 
we neglect intermittency and use the average dissipation 
rate, and then we include the effect of intermittency using 
the log-Poisson model. Assuming a v = 2 km s" 1 as an 
example, we find that the temperature increases to 63 K, 
much smaller than the 1000 K required for production of 
CH + . Temperatures much higher than the average can 
only be obtained if intermittency is included. In the tem- 
perature range from 100 K to 8000 K, the cooling rate 
as a function of temperature is very flat, meaning that 
hotter gas is more difficult to cool. This suggests that 
temperature fluctuations in this range of values reduce 
the overall cooling efficiency. Therefore, the net effect of 
fluctuations in the dissipation rate is that of producing 
a higher mean temperature than in the case of a uni- 
form dissipation rate (contrary to the case of molecular 
clouds, where the cooling rate depends on temperature 
very sensitively). For example, if the dimension of the 
most intense dissipative structures is d = 1.64, the av- 
erage temperature derived for er v = 2 km s" 1 is 100 K 
(versus 63 K in the case of a uniform dissipation rate). 
Furthermore, the flat cooling curve helps the generation 
of hot regions due to the fluctuations in the dissipation 
rate, so the extended high-rate tail of the dissipation rate 
distribution translates into a fairly extended high tem- 
perature tail in the temperature distribution. 

We calculated the gas temperature probability distri- 
bution using d — 1.64 for the most dissipative struc- 
tures in the intermittency model. Figure [7] shows the 
cumulative probability of temperatures lower than T for 
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different rms velocities. From the cumulative probabil- 
ity curve, we can read the fraction of gas with temper- 
ature larger than 1000 K. For a v = 1, 2, and 3 km s -1 , 
the probabilities of finding T > 1000 K are 3.6xl0~ 5 , 
3xl0 -3 and 0.019, respectively. Since the hot gas frac- 
tion needed to explain the observed abundance of CH + 
molecules is a few percent (Gredel et al. 1993), our re- 
sult shows that cold neutral clouds, where molecules such 
as CH + are detected, must have a minimum rms turbu- 
lent velocity of 3 km s _1 at the scale of 1 pc 3 . The 
uncertainty in d within the range 1.56-1.71 corresponds 
to a range of the cumulative probability of only ~ 2- 
3% around the value of 0.019 found for d = 1.64 and 
er v = 3 km s _1 . Therefore, the range of d does not affect 
our conclusion for the minimum rms velocity required to 
produce enough hot regions to explain the abundance of 
CH + molecules. One may expect the large uncertainty 
at the far right tail (Fig. 5) may significantly affect the 
cumulative probability. However, for a v = 3 km s" 1 , 
1000 K corresponds to a dissipation rate only 10 times 
larger than the average (see below for a discussion of a 
wider range of d, between 1 and 2 ). Therefore the cu- 
mulative probability is not determined by the far tail, 
but mainly by the part of the pdf where there are no sig- 
nificant snapshot-to- snapshot variations, as can be seen 
from Fig. [5] Figure [7] also shows that most of the gas 
has temperature below 100 K, which is the mean tem- 
perature in the case of a v — 2 km s _1 . This is because 
the skewness of the dissipation rate distribution favors 
values below the average. 

To study how the degree of intermittency affects the 
mass fraction of hot regions, we consider different val- 
ues (beyond the range 1.56-1.71) of the dimension of the 
most intense dissipative structures. We choose two spe- 
cial vales, d = 1 and d = 2, corresponding to vortex tubes 
and shocks. The dissipation rate is more intermittent for 
larger values of d. We find that, if a v = 1 km s _1 , the 
cumulative probability for T > 1000 K depends on d very 
sensitively It is 2 x 10~ 6 and 2 x 10~ 4 , respectively, for 
d = 1 and d — 2. However, for a v — 3 km s _1 , the 
degree of intermittency does not affect the probability 
by much; the cumulative probability is 0.015 and 0.021 
for d = 1 and d — 2, close to the value of 0.019 that 
we previously found for d — 1.64. This is because for 
<7 V = 1 km s _1 the average dissipation rate is low, and 
a temperature of T = 1000 K requires the far tail of 
the dissipation rate distribution, which strongly depends 
on the degree of intermittency. For er v = 3 km s _1 , in- 
stead, T = 1000 K corresponds to a dissipation rate only 
approximately 10 times larger than the average. In the 
range of values about an order of magnitude around the 
average, there is no considerable difference in the dissi- 
pation rate distribution for different dimensions of the 
dissipative structures. 

3 We note that Sheffer et al. (2008) find a similar minimum ve- 
locity dispersion from a different model proposed by Federman et 
al. (1996) to explain the observed CH+ abundance and its decou- 
pling from the OH abundance. In their model, reactions between 
neutral and ionic molecules, e.g., the production of CH+ from C+ 
and H2, are accelerated by the neutral-ion velocity, which is as- 
sumed to be the amplitude of MHD waves and is accounted for by 
using a non-thermal exponent for the effective temperature. We 
will investigate the intermittent distribution of the neutral-ion ve- 
locity in MHD turbulence and its effect on chemistry in HI clouds 
in a future work. 



Falgarone and Puget (1995) adopted an rms velocity 
of 3 km s _1 at 1 pc, and used experimental results for 
the probability distribution of the dissipation rate of in- 
compressible turbulence (converted from the velocity dif- 
ference distribution). Their model corresponds to our 
model with d — 1 . They obtained a mass fraction of a few 
percent (see their Figure 4) for gas with T > 1000 K, con- 
sistent with our result using the same parameters. There- 
fore, Puget and Falgarone's model, physically valid only 
if the turbulence in cold atomic clouds were subsonic, 
gives by chance a result similar to that of our model for 
supersonic turbulence. 

In summary, we have shown that the intermittent tur- 
bulent heating can generate enough hot regions to pro- 
duce the observed CH + abundance in cold neutral clouds 
if the turbulent rms velocity, cr v , is large enough. The re- 
quired minimum value for tr v is approximately 3 km s _1 
at a length scale of 1 pc. 

5. DISCUSSION AND CONCLUSIONS 

We have studied the energy dissipation and the heat- 
ing in supersonic turbulence. The turbulent dissipation 
is characterized by strongly intermittent fluctuations. A 
significant fraction of the kinetic energy is viscously dis- 
sipated in the finest, most intermittent structures, giving 
rise to a broad tail in the probability distribution func- 
tion of the dissipation rate. To study the turbulent heat- 
ing in interstellar clouds, a theoretical model is needed for 
the probability distribution of the dissipation rate at the 
dissipation scale, 77. We have generalized the log-Poisson 
model, originally proposed for incompressible turbulence 
by She and Leveque (1994), to supersonic turbulence. 
Because the dissipation scale, 77, cannot be resolved by 
current numerical simulations, we have used results from 
resolved inertial-range scales in our numerical simulation 
of supersonic and super- Alfvcnic turbulence as a guide- 
line for the sub-grid scales. We have found that the log- 
Poisson model, with a fractal dimension d — 1.64 for the 
most intermittent dissipative structures, gives an excel- 
lent fit to the mass-weighted probability distribution of 
the dissipation rate at resolved scales in the simulation. 
Extrapolating the model to the corresponding Reynolds 
numbers, we have studied the turbulent heating in molec- 
ular clouds, high-mass star-forming cores, and cold dif- 
fuse neutral clouds. Here we summarize our results. 

1. In typical molecular clouds, the average turbulent 
heating rate exceeds the cosmic ray heating rate 
by a factor of 3-4. Fluctuations in the heating rate 
give a lower mean temperature than expected from 
the average heating rate; temperature fluctuations 
make cooling more efficient because in molecular 
clouds the cooling rate increases sensitively with 
temperature. Taking intermittency into account, 
the turbulent heating alone gives a mean tempera- 
ture of approximately 8.5 K, close to the character- 
istic temperature of molecular clouds. This would 
suggest that cosmic rays are not even needed to 
explain the thermal balance in molecular clouds. 
However, due to the intermittent fluctuations in the 
turbulent heating rate, a significant mass fraction 
of the molecular gas is not heated by the turbu- 
lence. Cosmic ray heating would dominate in these 
regions. Assuming a cosmic ray heating rate of 
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0.8 x 1CP 2 n ergs cm -3 s -1 , the turbulent heating 
increases the average temperature by a few degrees, 
from 9.5 K to 13 K. We also find that turbulent 
heating plays a similar role in the outer regions of 
the dark cloud cores in the Orion molecular cloud 
complex. 

2. Turbulent heating provides an important energy 
source for the molecular gas in high-mass star- 
forming cores. Assuming spherical symmetry, in 
the inner regions of these cores (within ~ 0.01 pc 
from the center), stellar sources heat the dust 
grains to a relatively high temperature. Due to the 
high density of these inner regions (~ 10 6 cm -3 ), 
the gas and dust are thermally coupled, and the 
gas temperature is close to the dust temperature 
even in the absence of turbulent heating. The tur- 
bulent heating increases the gas temperature only 
by a few K. On the other hand, in the outer re- 
gions (~ 0.1 pc from the center), where most of the 
core mass resides, the turbulent heating results in a 
considerable increase in the gas temperature. The 
low density of these regions (10 4 -10 5 cm -3 ) makes 
energy exchange between gas and dusts inefficient 
and, without turbulent heating, the energy transfer 
from the dust heats the gas only to 10-20 K. Inclu- 
sion of turbulent heating increases the temperature 
to approximately 36 K. Because turbulent heating 
causes a large temperature increase in most of the 
core mass, it may have important implications for 
the dynamical evolution of the cores and for their 
star-formation process, and it may also be probed 
by future observations. 

3. The intermittent turbulent heating in diffuse HI 
clouds can give rise to regions much hotter than 
the average temperature. These warm regions have 
been used to explain the existence of molecules such 
as CH + in HI clouds, whose production needs high 
temperatures. Using the log-Poisson intermittency 
model for supersonic turbulence, we find that a 
turbulent rms velocity of 3 km s _1 at 1 pc is suf- 
ficient to account for the observed abundance of 
these molecules, which extends the earlier result, 
based on incompressible turbulence, by Falgarone 
and Puget (1995). 

We point out that thermal conduction is neglected in 
our calculations. Conduction tends to transport ther- 
mal energy to fill in regions not significantly heated by 
dissipation, and thus may, to some degree, erase the fluc- 
tuations in the heating rate. As discussed in § 4, more in- 
termittent fluctuations in the heating rate give a smaller 
average temperature in molecular clouds. Therefore, if 
thermal conduction were included, the average temper- 
ature in molecular clouds and in the outer regions of 
high-mass star-forming cores would be even larger, mak- 
ing turbulent heating even more important in these two 
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cases (§ 4.1 and § 4.2). The situation is different in cold 
diffuse HI clouds (§ 4.3). If the fluctuations in the heat- 
ing rate are less intermittent there, the average temper- 
ature becomes smaller and, more importantly, the tail 
of the temperature probability distribution, needed for 
CH+ production, would be less extended. We estimate 
whether and how much thermal conduction would change 
our results by calculating the conduction length scale, 
l c ~ \fnt c , during a cooling time scale, t c , where k is the 
thermal conduction coefficient (approximately equal to 
the kinematic viscosity, v). This is the scale over which 
thermal conduction can homogenize before the heat from 
turbulent dissipation is radiated away. We find that in 
molecular clouds and in HI clouds l c is smaller than (but 
comparable to) the dissipation length scale, rj. There- 
fore, the heat generated in the most intermittent struc- 
tures cannot be transported far from these structures by 
thermal conduction. This justifies our choice of comput- 
ing the temperature distribution by using the dissipation 
rate distribution evaluated at r\. In the outer regions of 
high-mass star-forming cores, l c is about an order of mag- 
nitude larger than rj (the latter is very small because of 
the very large Reynolds number). Using the distribu- 
tion of the dissipation rate at l c instead of rj (assuming 
thermal conduction homogenizes the temperature over a 
size of l c ) gives an average temperature of 40 K in the 
outer regions of these cores, a little higher than from the 
distribution of the dissipation rate evaluated at 77. How- 
ever, the conduction process in the presence of turbulent 
motions is more complex than described by the above 
estimate, and a detailed study of its effects is beyond the 
scope of this paper. 

Although so far neglected, the process of turbulent 
heating may play an important role in the process of star 
formation. The mean temperature in molecular clouds 
defines the mean Jeans mass, which may control the 
peak of the stellar mass distribution. Because we have 
found that the mean temperature in molecular clouds 
and in high-mass star-forming cores may be controlled 
by turbulent heating, the characteristic stellar mass may 
be affected by turbulent heating as well. For example, 
the larger gas temperature predicted in high-mass star- 
forming cores may partly offset their large density, result- 
ing in almost the same characteristic stellar mass, with 
respect to molecular clouds following Larson's relations. 
Finally, because the process of turbulent heating results 
in broad gas temperature distributions, it may be cru- 
cial in many molecular chemical reactions, besides those 
responsible for the formation of CH + molecules. 
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